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PREFACE 

Precipitation  static  has  long  been  recognized  as  a  particularly  serious 

problem  in  aviation  radio  navigation  and  co.amunication  systems.  Many  investi- 

8,  11,  20 
gations  as  to  its  causes  and  nature  have  been  made        with  successful 

progress,  especially  during  the  past  fifteen  years.  Static  free  navigational 
systems  (VDR  Omni-range),  for  short  ranges,  and  line  of  sight  corfimunications 
systems  (VHF-UHF),  are  now  standard  equipment  on  all  aircraft,  although  it 
must  be  pointed  out  that  these  solutions  involve  a  sort  of  strategic  re- 
treat, that  is,  a  shift  to  higher  frequency  levels  where  the  precipitation 
static  energy  content  is  of  a  lesser  value.  Also  the  investigations  that 
were  made  all  utilized  propeller  driven  aircraft  for  their  test  vehicles 
and  although  projects  were  earmarked  to  include  jet  type  of  aircraft,  none, 
to  the  author's  knowledge,  were  ever  reported  on. 

The  problem  that  is  the  subject  of  this  thesis  was  essentially  to  de- 
sign a  device  to  measure  the  relative  noise  power  or  voltage  due  to  corona 
discharge,  at  certain  specified  locations  in  the  frequency  spectrum,  and 
be  carried  aboard  jet  type  aircraft.  This  is  a  bit  of  a  swift  explanation 
but  a  detailed  account  will  be  given  later. 

Also  an  attempt  is  made  to  correlate  information  regarding  some  of 
the  statistical  aspects  of  noise  analysis  and  although  it  is  not  a  completely 
rigid  treatment,  it  describes  in  a  condensed  manner,  a  description  of  noise 
signals  and  their  effects  on  electronics  circuits. 

The  author  would  like  to  express  great  appreciation  to  Messrs.  H.  P. 
Blanc hard  and  Maurice  Mills  of  the  Stanford  Research  Institute  for 
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their  infinite  patience  and  guidance,  and  to  Professors  Robert  Bauer  and 
Robert  Kahal  for  encouragement  and  interest. 
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Chapter  I 
INTRODUCTION 

1.  History  and  General  Considerations  of  the  Problem 

Precipitation  static  may  be  described  as  radio  interference  caused  by 
corona  discharge.  It  is  sometimes  labled  as,  'Poisson  noise',  is  tech- 
nically named,  'impulse  noise',  and  is  statistically  described  as,  'nearly 
normal  random  noise'.  It  is  but  a  single  type  of  the  many  kinds  of  noise 
that  have  been  studied,  analyzed  and  catalogued,  although  it  does  have  the 
dubious  distinction  of  being  the  number  one  culprit  of  aviation  naviga- 
tional and  communications  systems. 

The  first  analysis  and  explanation  of  precipitation  static  was  that 
charged  particles,  created  by  the  turbulence  and  charged  energy  of  thunder 
clouds,  generated  static  when  striking  the  antenna.  No  attempts  were  made 
to  alleviate  these  effects  for  many  years  until  the  introduction  of  the 
shielded  loop  antenna  in  193U*  It  was  assumed  that  since  particles  did  not 
strike  the  antenna  proper,  static  would  be  reduced.  However,  this  assump- 
tion was  not  substantiated  in  practice  because  it  was  found  that  in  actual 
air  operation,  severe  static  easily  disabled  the  system. 

In  1939  after  more  investigation,  a  new  theory  ^  assigned  the  main 
cause  of  precipitation  static  in  aircraft  to  that  of  corona  discharge. 
Later  the  concentrated  effort  of  the  Army-Navy  Precipitation  Static  Project  ° 
verified  this  theory  which  resulted  in  the  subsequent  changes  to  air-frame 
construction  and  airborne  electronic  equipment  that  are  in  standard  prac- 
tice today. 

With  the  increased  use  of  jet  and  rocket  propulsion  in  supersonic 
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aircraft  and  guided  missiles,  it  was  erroneously  assumed  that  the  determi- 
nental  effects  due  to  precipitation  static  would  be  of  no  consequence, 
since  in  this  type  of  craft,  a  natural  discharging  process  due  to  the  ther- 
mal ionization  of  the  exhaust  gases,  should  efficiently  transfer  free  charge 
to  the  atmosphere.  This  is  partially  accomplished  because  a  complete  recom- 
bination of  the  ions  produced  in  burning  the  fuel  has  not  as  yet  taken  place, 

thus  the  electric  field  of  the  aircraft  will  capture  those  ions  of  opposite 

p 
sign  while  simultaneously  repelling  those  of  the  same  sign. 

In  a  recent  report  by  R.  L.  Tanner,  °it  was  shown  that  the  strength 
of  the  noise  signals  due  to  corona  discharge  is  directly  related  to  the 
dipole  moment  of  the  discharge  current.  The  dipole  moment  is  proportional 
to  the  product  of  the  current  flow  and  the  distance  over  which  it  takes 
place.  Thus  in  any  supersonic  aircraft  which  strikes  ice  crystals  (cumulo- 
cirrus  type  clouds),  or  heavy  turbulence  (cumulo-nimbus  type  clouds),  or 
the  layers  of  ionized  air  located  at  extremely  high  altitudes,  a  static  charge 
builds  up  in  the  vicinity  of  plastic  antenna  housings,  canopies  and  control 
surfaces.  It  is  here  that  one  is  faced  with  a  new  type  phenomenon  that  is 
based  on  a  very  old  principles  that  substances  of  high  dielectric  constant 
acquire  a  positive  charge  when  rubbed  with  substances  of  lower  dielectric 
constant.  This  is  called  triboelectric  charging,  and  can  be  visualized  as 
a  scrapping  off  of  the  electrons  by  the  air stream.  Naturally  this  charging 
eventually  culminates  in  a  corona  type  discharge  at  many  points  close  to 
the  outer  perimeter  of  plastic  housings  and  the  dipole  moment  of  the  cur- 
rent discharge  acts  as  an  efficient  antenna  to  couple  this  electrical 
noise  into  the  receiver  input  circuits. 

The  problem  specifically,  is  to  determine  the  average  noise  power  or 
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the  average  voltage  produced  by  this  static,  recorded  for  a  period  of 
time  equal  to  the  flight  duration  of  the  aircraft.  This  system  must  dupli- 
cate an  actual  navigational  system's  characteristics  with  respect  to  band- 
width, resonant  frequency,  servo  system  time  constants  and  antenna  match- 
ing stage. 

The  actual  navigational  system  can  be  represented  in  general  by  the 
diagram  in  Figure  1. 
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FIGURE  1 
Simplified  Diagram  of  Navigational  System 


To  illustrate: 

If  the  incoming  waveform  is  sinusoidal  in  nature,  then  the  output  of 
the  mixture  will  be  the  sum  and  difference  terms  of  the  original  inputs. 
This  is  applied  to  the  integrator  which  acts  to  average  only  the  difference 
terras  since  it  is  a  low  pass  filter.  The  output  will  then  be  a  very  slowjLy 
varying  voltage  and  the  subsequent  error  signals  will  also  be  very  small 
in  magnitude. 

Where  precipitation  static  is  present  the  waveforms  (idealized)  might 
take  the  form  of  those  shown  in  Figure  2.  Since  power  is  proportional  to 
the  voltage  applied  to  the  servo  system,  the  disturbance  to  the  system  is 
related  to  all  the  input  noise  pulses  integrated  over  a  period  of  time 
equal  to  the  time  constant  of  the  system,  and  can  be  written  as 


— A  C-  fo%  —  Phase  Disturbance  in  System 


The  error  signal  to  the  reinsertion  oscillator  will  then  be  proportional 
to  the  input  average  noise  power  so  that  periods  of  heavy  static  would  pre- 
cede large  error  signals. 

The  system  that  was  designed  and  constructed  duplicating  the  above 
characteristics  is  now  being  operationally  used  to  record  the  average  noise 
power  over  periods  of  ten  to  twelve  hours.  These  results  are  later  com- 
pared to  the  navigational  errors  plotted  separately  by  the  navigator  to  see 
what  relative  values  of  average  noise  power  cause  the  maximum  errors  in 
the  system. 
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FIGURE  2 
Idealized  Waveforms 
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FIGURE  2 
Idealized  Waveforms 


A  second  system  to  measure  the  average  noise  power  at  specific  fre- 
quencies is  also  described  in  Chapter  III  with  some  data  which  shows  re- 
sults of  actual  flights  made  using  the  system. 


2.  The  Projected  Method  of  Solution: 

A  study  of  the  available  systems  for  measuring  noise  due  to  atmospheric 
sources  and  corona  discharge  indicates  a  common  tendency  to  make  the  major- 
ity of  measurements  using  an  electric  field  intensity  meter,  or  circuits  of 
some  such  similar  nature.  The  results  are  then  interpreted  using  a  graph- 
ical portrayal  of  corona  discharge  current  as  a  function  of  electrical  field 
intensity.  In  this  problem  it  is  desired  to  measure  the  energy  components 
of  the  noise  at  a  specified  frequency  and  known  bandwidth,  therefore  a 
relationship  with  respect  to  time  must  be  set  up  to  measure  this  noise  power. 
But  the  instantaneous  noise  voltage  induced  in  the  antenna  cannot  be  con- 
sidered as  a  regular  function  of  time,  for  if  it  were,  both  amplitude  and 
phase  relationships  would  have  to  be  known  for  one  complete  period.  Since 
it  cannot,  the  problem  is  one  that  is  statistical  in  nature. 

Statistics  may  be  succinctly  defined  as  the  science  of  reduction  and 
analysis  of  data.  Its  application  to  the  problems  associated  with  com- 
munication processes  can  be  mainly  attributed  to  N .  Wienerf ^  Who  was  one 
of  the  foremost  contributors  during  the  development  years  of  statistical 
prediction  theory,  particularly  concerning  stationary  random  time  func- 
tions, a  random  time  function  being  defined  as  a  function  which  does  not 
follow  any  specific  law  of  variation  and  therefore  is  not  precisely  predic- 
table. The  type  of  noise  being  investigated  in  this  paper  falls  into  this 
category  for  it  consists  of  a  superposition  of  impulses  occurring  at  ran- 
dom times  with  overlapping  individual  pulses.  It  is  sometime  referred  to 
as  Poisson  noise  because  the  Poisson  distribution  (Appendix  I)  is  used  to 
describe  its  statistical  structure.  In  statistical  applications  the 
Poisson  Distribution  often  appears  when  concerned  with  the  number  of 
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occurrences  of  a  certain  event  for  a  very  large  number  of  observatione, 

3,5 
where  probability  for  the  event  to  occur  in  each  observation  is  very  small. 

It  is  for  this  reason  that  it  is  sometimes  called  the  "small  probability  law". 

For  the  case  of  very  high  density  Poisson  noise,  that  is,  when  the  number 

of  elementary  random  impulses  is  large,  and  their  mean  duration  is  not 

too  brief,  a  heavy  overlapping  occurs  and  the  result  is  nearly  normal  ran- 

17 
dora  noise.  This,  according  to  Middleton  '  can  be  considered  to  have  a 

normal  distribution  with  negligible  correction  terms,  and  the  shape  of  the 

pulse  is  unimportant.  (Appendix  I) 

It  is  plausible  to  assume  that  the  corona  discharge  that  takes  place 
will  have  the  characteristics  of  high  density  Poisson  noise,  for  when  the 
same  navigational  system  was  operated  in  slower  propeller  driven  types  of 
aircraft,  the  system  suffered  erratic  operation  (intolerable  errors), 
only  when  thunder  storms  set  up  exceedingly  heavy  corona  conditions.  Also, 
the  system  operated  normally  under  conditions  of  moderate  corona  indicating 
that  the  low  density  types  can  be  ignored. 

For  the  case  of  the  jet  aircraft,  the  corona  discharge  intensity  would 
be  mainly  a  function  of  speed  and  since  the  latter  is  usually  close  to 
maximum  designed  speed,  the  intensity  will  always  be  close  to  maximum, 
with  corona  discharge  pulses  occurring  at  all  of  the  places  where  plastic 
surfaces  are  exposed  to  the  airstream.  The  mathematical  theory  to  des- 
cribe this  phenomena  is  covered  below  and  is  closely  allied  to  the  more 
familiar  treatment  of  fluctuation  noise. 


3»  A  development  of  Mathematical  Description  of  Random  Noise. 

Consider  a  plane-electrode  diode  for  the  emission  limited  case,  A 
graph  showing  the  current  versus  time  for  a  single  electron  leaving  the 
point  of  zero  potential  (  t  =  0  ),  and  arriving  at  the  point  of  highest 
potential  (  t  =  tn),  is  shown  in  Figure  3.  When  a  constant  potential  is 
assumed  over  a  time  interval  2T,  and  a  pulse  having  the  shape  of  that 
shown  in  Figure  3>  takes  place  for  each  electron  crossing,  a  plot  of 
the  current  l(t)  versus  time  for  the  interval  2T  is  shown  in  Figure  1*. 

Assume  that  N  electrons  are  emitted  during  the  interval  indicated. 
Then  l(t)  is  the  sum  of  the  individual  current  pulses 
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and  the  time  average  is 
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FIGURE  3 
Graph  of  i(t)  versus  Time 
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FIGURE  4 
Graph  of  i(t)  versus  Time  Showing  the  Interval  2T 
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Now  consider  the  electron  stream  to  be  made  up  of  a  great  number  of 

19 
electrons  whose  fluctuations  are  considered  to  be  random.    A  count  of 

electrons  arriving  at  the  plate  in  the  interval  2T  is  N.  Repeating  this 

count  for  a  total  number  of  M  times  gives  a  value  for  the  average  number 

of  electrons  per  second  as 


As  M  is  increased  and  T  held  constant,  some  of  the  N's  will  have 
the  same  value,  thus  the  number  of  N's  of  a  particular  value  will  in- 
crease as  M  is  further  increased.  The   probability  of  getting  N  electrons 
for  some  particular  trial  is  defined  as 

p(N)  =  lim   Number  of  trials  giving  N  electrons  \j$ 


Since  p(N)  also  depends  on  T  it  can  be  assumed  that  the  randomness  of 
the  electron  stream  is  such  that  the  probability  that  an  electron  will 
arrive  at  the  anode  in  the  interval  (t,   t  ■*•  A  t)  is 
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v^t  where  At  is  such  that  vAt-*<l,  and  that  this  probability  is  in- 
dependent of  what  happened  before  time  t  or  what  will  happen  after  time 
t+At.     For  a  very  large  number  of  trials  this  is  a  description  of  the 
Poisson  distribution  and  is  expressed  as 

p(N)     «  (2W)N  e"2vT  X-6 

K 

Tnus  a  large  number,  M,  of  intervals  of  length  2T  having  exactly 
N  arrivals  will  be  to  a  first  approximation  Mp(N).  For  a  fixed  value  of 
t,  equation  1-3  may  be  written 

and  for  A<"t^  T~"&     it  effectively  becomes 

dt 

1-8 


5»  -  %\  to) 

\J-oc 


It  is  now  possible  to  write  the  average  l(t)  over  all  of  the  M  in- 
tervals instead  of  only  over  those  having  N  arrivals,  so  as  M 


oo 


so 
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Campbell's  theorem1^  states  that  the  mean  square  value  of  the  fluctuation 
about  the  average  is 

,00 


.. 


(I(t)  -  I(t)T  =       ▼     \     i  U)  dt  1-12 

•o* 


i2u)  -  2i(t)  TTT5      -   i(t)*         1-13 


Since 


(t)    s  £  £     i(t  -  tft)  i(t  -  tm)  i-iii 

A*\  rr\ts\ 

Averaged  over  all  values  of  t^>  ^>2>  ^3>      ^N   »  with  t  held 

fixed 


The  multiple  integral  has  two  different  values,  so  tfien  n  ■  n 


zr 

and  for  n  ^  m  its  value  is 
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zr 
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IT  1-17 


Counting  up  the  number  of  terms  in  the  double  sum  shows  that  there  are 
N  of  them  having  the  first  value  and  N2  -  N  having  the  second  value. 


1U 


For       -t  +  ^<   t<^     T  -  A 

1-18 


i>>-fr&^^rr,-«* 


Averaging  over  all  the  intervals  instead  of  only  those  having  N 
arrivals,  the  final  form  can  be  written  as 


x&  - 1  ^  i;cf) 


1-19 
KJ=0 

poo  

1-20 


where  the  summation  with  respect  to  N  is  performed  as  in  (l-ll)  and  after 
the  summation,  the  value  (1-11)  for   I(t) is  used. 

The  relationship  of  the  above  development  to  the  problem  can  be  shown 
for  a  small  increment  of  surface  area  located  at  one  of  the  points  where 
the  metallic  airframe  and  plastic  housing  are  adjacent. 

By  equation  1-1  the  individual  current  pulses  can  be  described  for  a 
single  element,  and  when  all  of  the  points  are  taken  into  consideration, 
plus  all  of  the  locations  where  plastic  housings  are  adjacent  to  metallic 
airframe,  equation  1-20  expresses  the  total  average  squared  current. 
(Appendix  I) 
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1*.  Previous  Statistical  Noise  Measurements 

Some  investigation  of  other  statistical  approaches  to  the  measure- 
ment of  noise  was  made  before  any  design-  considerations  were  conceived. 

In  their  article,  "A  Statistical  Approach  to  the  Measurement  of  At- 
mospheric Noise"  by  R.  S.  Hoff  and  R.  C.  Johnson,   it  states  in  part,"... 
that  atmospheric  noise,  causing  noise  voltage  to  be  induced  into  a  receiv- 
ing antenna.. .the  envelope  of  this  noise  voltage  after  being  amplified 
by  the  receiver  stages  and  further  modified  by  limiting  and  other  non- 
linear processes,  is  an  irregular  function  of  time  or  time  series  I(t), 

Ik 
and  may  be  treated  by  statistical  processes."  Also,  Landon   has  shown 

that  the  distribution  of  noise  voltage  when  passed  through  a  narrow  band 

filter  is  normal  (Gaussian), (Appendix  I)  that  is,  the  fluctuation  noise 

output  of  a  band  pass  amplifier  has  a  distribution  of  voltage  versus  time 

which  follows  the  normal  error  law.  The  method  used  is  similar  to  that  used 

-I  o 

by  Lord  Rayleigh   in  determining  the  probable  distribution  with  time  of 

the  instantaneous  amplitude  of  a  large  number  of  sine  waves,  all  with 

12 
equal  amplitudes  and  random  relative  phases.  From  Knudzton's  report  ,  "... 

that  the  narrow  band  noise  has  the  character  of  an  amplitude  and  frequency 
modulated  wave  for  the  ideal  rectangular  band  pass  filter,  where  the 
average  period  of  the  envelope  is  a  linear  function  of  the  bandwidth  only." 

Thus  it  is  possible  to  list  the  conditions  that  will  effect  the  de- 
sign characteristics  of  the  device  to  examine  the  properties  of  the  pre- 
cipitation noise, 

(1)  If  a  quantitative  value  of  the  noise  power  is  desired,  a  narrow 
band  pass  filter  of  approximately  3  kc  must  be  used  since  the  center 
frequency  is  126  kc .  (See  Appendix  II) 

(2)  The  device  must  simulate  a  navigational  system  which  has  the  fol- 
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lowing  characteristics: 

(a)  Bandwidth  30  kc 

(b)  Resonant  Frequency  126  kc 

(c)  Servo  time  constant  response  -  100  seconds 

(3)  Ideally  considered,  the  noise  has  a  white  spectrum.     This  of  course 
is  not  true,  but  is  permissible  for  this  problem  since  Plank's  radiation 
formula  indicates  that  the  thermal  noise  energy  is  down  l/e  of  its  maxi- 
mum value  of  KT  at  about  10'  mc.^ 

To  satisfactorily  solve  this  dilemma  of  having  two  bandwidths  it  would 
be  necessary  to  have  two  distinct  systems;  the  first  to  measure  a  quanti- 
tative value  of  the  noise,  and  the  second  to  ride  next  to  the  navigational 
system  recording  the  relative  noise  input  level  for  a  specified  naviga- 
tional track.     It  was  decided  to  use  the  latter  system  to  study  the  effects 
peak  disturbance  periods  as  a  function  of  navigational  error  and  to  use 
the  first  system  to  set  up  a  method  to  investigate  the  power  spectrum  of 
high  density  precipitation  static  by  measuring  the  average  power  level 
across  the  band  in  small  increments  of  frequency. 
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Chapter  II 
Description  of  Designed  Measuring  Devices 

1.   SYSTEM  I  -  MEAN   SQUARE  NOISE  METER 

As  shown  in  Appendix  II,   it  is  possible  to  represent  the  output  of 
a  narrow  band  filter  by 

e(t)     s     e  (t)  cos  w„t  4*   e   ,      (t)   sin  w  t 

cos  m     ~      sin  v    '       "    m 

Using  a  very  narrow  band  pass  amplifier  as  the  input  stage,  the  mean 
value  of  the  output  noise  voltage  can  be  found  by  integrating  the  voltage 
over  a  long  period  of  time.  This  integrator  has  an  output  wiach  is  made 
up  of  a  component  that  increases  linearly  with  time,  plus  some  variation 
due  to  the  input  a-c  components  wnich  are  smoothed  out  after  sufficient 
elapsed  time.  It  is  specifically  desired  to  record  the  mean  square  value 
of  the  noise  voltage  (  e*),  by  constructing  a  device  which  will  square  the 
input  and  then  integrate  the  low  frequency  (difference  terms)  components 
over  a  long  period  of  time.  The  advantages  are: 

(a)  The  recorded  measurements  are  proportional  to  energy. 

(b)  From  the  basic  premise,  that  the  noise  has  a  normal  distribution, 
it  must  then  consist  of  a  uniform  spectrum  of  Spectral  Density  N 

volts  /cps  at  all  frequencies  other  than  zero,  and  have  a  true  mean  value 
(e)volts,  it  is  necessary  to  continue  the  integration  process  for  a  length 
of  time 

T  =  (2  x  lcA  NQ)  /P2  e*  2sec 
in  order  to  insure  a  95%   certainty  that  the  measured  value  lies  within 
P  percent  of  the  true  mean. 
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Frequency  Analysis  of  the  Square  Law  Circuit:     The  current  of  a 

diode  may  be  represented  by  the  infinite  series 

2            3  n 

ip  -  t^e  +•  c2  e  ■*  cyrV *cne 

in  which  the  constants  cn  may  be  determined  as  are  the  coefficients  in  a 
Taylor  series. 

Since  the  diodes  are  operated  in  a  region  where  the  dynamic  character- 
istic is  not  a  straight  line,  let  it  be  assumed  that  the  first  two  terms 
accurately  represent  the  portion  of  the  characteristic  used.  Thus 

S  "   cl  e+c2  e* 

Substituting  equation  III-7:  Appendix  II 

S>  =  Cl  (  ec^   cos  V  +  esU)  sin  wmt) 

+  *2   (   ec^   cos  wmt'+es(t)  sin  vml   )2 

The  only  terms  of  interest  are  those  containing  difference  terms  since 
the  integrator  is  essentially  a  low  pass  filter,  therefore  in  this  case 
it  is 

2  c2  ec(t)  eg(t)  sin  wm  t  cos  wmt 
By  trigonometric  identities  trds  is  equal  to 

2  c2  ec(t)  eg(t)  I  sin  (wm  -  wQ)  t  +  cos  wmt  sin  w0tj 
Thus  the  presence  of  the  difference  terms  allows  the  use  of  a  large 
time  constant  in  the  integrator  circuit  to  average  the  input  low  fre- 
quency terms. 
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Description  of  Circuits  in  System  I 

Figure  1-0     Block  Diagram  of  the  System 

Figure  1-1    Narrow  Band  Pass  Amplifier  -  Essentially  a  tuned  grid, 
tuned  plate  tank,  measured  to  be  3kc  wide  at  half  power  points,  with  center 
frequency  at  126  kc,  and  tunable  over  a  30  kc  band. 

Figure  1-2  -  Logarithmic  Amplifier  -  The  amplifier  is  an  instan- 
taneous linear  amplifier  designed  to  obtain  a  large  linear  range.     A  feed- 
back loop  from  the  Miller  integrator  circuit  which  has  a  time  constant  of 
100  seconds  causes  a  compressing  effect  to  the  input  signal.     This  allows 
examination  of  signals  whose  magnitude  may  vary  from  10  micro-volts  to 
10^  micro-volts.     The  variable  gain  is  calibrated  to  give  approximately 
10  db  gain  for  -  6  volts  of  avc  and  70  db  gain  for  -  1  volt  avc. 

Figure  1-3     Square  Law  Circuit  -  This  circuit  was  originally  designed 
using  the  relation 

E0  =  k  E±2 

where  EQ  is  the  output  voltage  and  E?  is  the  input  voltage  of  the  circuit. 
Knowing  the  desired  avc  characteristic  and  the  time  constant  of  the  Miller 
integrator  it  is  possible  to  approximate  the  value  of  k.  This  was  accom- 
plished by  plotting  several  curves  using  values  close  to  the  approximated 
value  of  k,  constructing  a  circuit  with  variable  resistors  and  adjustable 
bias  levels,  empirically  determining  the  component  values  after  the  input 
and  output  values    'fit'   the  desired  curve. 

Figure  1-k     Detector  and  Miller  Integrator  -  The  detector  is  a  peak 
reading  device  with  a  time  constant  determined  by  the  center  frequency 

T  C  detector    =     1/1.26  x  10*  =  8  micro-seconds 
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The  basic  principle  of  the  Miller  integrator  in  this  system  is  the 
operation  of  averaging.  It  is  performed  by  passing  the  general  input  through 
a  low  pass  filter  (integrator),  thus  the  contributions  to  the  output  will 
be  made  only  by  those  terms  which  have  frequencies  which  fall  within  the 
range  of  the  filter.  In  this  problem  the  observation  interval  is  fixed 
at  100  seconds,  so  thus  the  output,  M(T),  of  the  averaging  filter  may  be 
related  to  its  input  by  means  of  the  convolution  integral  involving  the 
filter's  impulse  response  and  weighting  function,  h(t), 


Si 


T 
M(T)     =         \  h(t)         E     (T  -  t)  dt 


where  T  is  the  duration  of  the  observation  interval. 

From  "Statistical  Errors  in  Measurements"   ,  Table  Ij   the  effective 
bandwidth,  Bjj  and  the  weighting  function,  h(t),  are  given  for  an  integrator 
as 

h(t)     =     l/T     =     o.Ol  where     0    >  t>T 

%         =     1/2T  s     0.005 
The  input  will  be  the  difference  terra 

Ed(t)     =     2  c2  ec(t)  es(t)  sin  (wm  -  wQ)t 

for 

0  <   w  <   0.01 

Simplifying  M(T) 
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The  above  represents  the  output  of  the  averaging  filter  for  an 
observation  interval  of  100  seconds. 

Figure  1-5  Meter  Coupling  Circuit  -  This  circuit  uses  an  Esterline 
Angus  type  recorder,  Model  AW.  The  meter  movement  is  a  1  ma  type. 
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2.  System  II  -  NAVIGATIONAL  DEVICE  SIMULATOR  &  RECORDER 

As  was  previously  stated,  this  system  must  duplicate  a  navigational 
device  that  is  represented  by  the  diagram  in  Figure  1.  It  must  record 
noise  voltage  but  must  disregard  the  signal,  so  that  if  the  circuit  shown 
in  Figure  1  were  used  in  its  entirety,  a  recording  of  the  output  would 
simply  be  the  average  value  of  the  input  signal  and  noise  lumped  together 
with  no  means  of  distinguishing  one  from  the  other.  It  is  required  then, 
that  the  signal  must  be  eliminated  and  only  the  noise  voltage  recorded. 
Filtering  of  course  will  be  useless  since  it  removes  those  noise  frequencies 
that  are  of  interest,  so  that  the  only  solution  is  to  phase  out  the  signal. 
This  will  require  two  separate  channels  and  the  use  of  phase  shifters,  but 
in  keeping  with  the  basic  premise,  to  duplicate  the  navigational  system 
as  much  as  possible,  the  only  change  of  significance  wiil  be  the  intro- 
duction of  a  balanced  modulator.  Its  purpose  is  to  ict  as  a  gate,  opening 
and  closing  at  a  set  rate,  but  90°  out  of  phase  with  the  gate  in  the  other 
channel.  Thus  when  the  two  channels  are  added  together,  the  signal  will 
have  been  phased  out  but  the  noise  components  will  be  the  sum  of  all  the 
cosine  terms  (from  one  channel)  and  all  of  the  sine  terms  (from  the  other 
channel) . 
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Description  of  Circuits  in  System  II 
Figure  II-O    Block  Diagram  of  System 

Figure  II-l    Band  Pass  Filter  and  Tuned  Amplifier  -  Essentially  tnis 
circuit  is  a  single  tuned,  link  coupled  circuit,  acting  as  a  band  pass 
amplifier  with  bandwidth  of  30  kc  with  a  resonant  frequency  of  126  kc. 
Figure  II-2   Logarithmic  Amplifier  -  Same  as  System  I. 
Figure  II-3   126  kc  Crystal  Oscillator  and  Phasing  Networks.  A  highly 
stable  Pierce  crystal  controlled  oscillator  is  used  as  a  local  oscillator, 
the  output  of  which  is  fed  to  two  RC  type  phasing  networks.  If  the  out- 
put of  the  oscillator  is 

eQ  ■  E  sin  w0t 
then  the  outputs  of  the  phasing  networks  will  be  respectively 

e-j^  «  E  sin  wQt  -  h$° 


e?  =  E  sin  w_t  +  1£° 


and 

^2  ~  ^  sin  wQi 

Since  e^  and  e.  are  ninety  degrees  apart  and  bear  no  fixed  relationship  to 

the  input  signal  as  yet,  they  can  be  represented  as 

e-,  =  E  sin  w  t 
J-        o 

&2  ~   E  cos  w0t 
Each  are  amplified  in  separate  channels  and  remain  in  these  channels  un- 
til detections  take  place. 

Figure  II-U   126  kc  Amplifier,  Gate  and  Integrator  -  The  first  two  stages 
are  voltage  gain  stages  which  split  e-^  into  two  voltages  180°  apart  in 
phase  and  thus  drive  a  push  pull  cathode  follower  stage.  By  this  means 
it  is  possible  to  offer  a  low  impedance  driving  source  to  the  balanced 
modulator  ( gate ) . 
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Frequency  Analysis  of  Coherent  Detector   (Balanced  Modulator) 
If  the  input  noise  voltage  has  the  form  at  the  output  of  the  log 
amplifier  of 

e(t)  s  Ax  (t)  sin  wt  +    A2   (t)  cos  wt 
where 

w  =  2  1T  ( 126  kc    *  30  kc ) 
the  output  of  the  balanced  modulator  can  be  computed  in  the  following 
manner. 

Assuming  matched  diodes  with  similar  current  characteristics,   the 
current  can  be  represented  by  the  infinite  series 

V     a0+aiera2e2t    ^     e^  +.  ....  +  an  en 

Since  the  integrator  will  act  as  a  low  pass  filter  on  the  difference 
terms  of  the  output  of  the  non-linear  device  will  be  of  consequence. 
The  d-c  term  (aQ)  can  be  ignored  since  transformer  coupling  is  used. 

The  voltage  at  point   (1),  Figure  II-1?  is 

%  a  e^eCt) 
The  voltage  at  point  (2),   Figure  II-£  ie 
Eg  s    .  e1^e(t) 

thus  the  value  of  e         can  be  written  as 

out 

=  i0  =  a2  e2  =    A(E2  -  e0)2  -  A(e0  -  E-J2 


R 
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Setting  AR  =  1/K 

K  eo  =     (-ei+  e(t)  -  eo)2  -  (eQ  -  ex  -  e(t)): 

Thus 

-1]^  e(t) 


K  -  h   ex 
Substituting 


eo  = 


-  liEAi  (t)  sin  w  t  sin  wt  -  h   E  A0(t)  sin  w  t  cos  wt 

-L_ O c.        O 

K  -  U  E  sin  w  t 


Since  only  the  difference  terms  are  of  consequence 

-  sin  v  t  cos  wt  =  -  sin  (w^  -  w)  t  -  cos  w^t  sin  wt 
o  O  0 

and 

sin  wQt  sin  wt  =  -  cos  (wQ  -  w)  t  -*-  cos  wQ  cos  wt 

Thus  eQ  can  be  represented  as  a  voltage  having  frequencies  made  up  of 

the  difference  terms. 

For  channel  two  the  analysis  is  the  same  with  the  exception  of  the 

value  of  e« 

-U  e2  e(t) 


eo  = 


K  -  U  e2 
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Substituting 

e  -  -U  E  A-l  (t)  cos  w0t  sin  wt  -  k  E  Ap(t)  cos  wQt  cos  wt 

K  -  k   E  cos  w  t 
o 

Again  only  the  difference  terms  are  of  consequence 

-  cos  vQt  sin  wt  «  sin  (wo  -  v)t  -  sin  w  t  cos  wt 

and 

-  cos  wQt  cos  wt  ■  gin  wQt  sin  wt  -  cos  (wQ  -  w)  t 

as  before,  the  voltage  applied  to  the  integrator  can  be  represented  as 
having  frequencies  made  up  of  the  difference  terms. 

The  Miller  integrator  is  similar  to  that  used  in  System  I  except 
that  the  output  can  be  either  a  positive  or  negative  value,  since  the 
average  will  vary  about  zero  and  will  have  a  mean  average  level  of  zero. 
Figure  II-£  U00  Cycle  Oscillator  and  Amplifiers  -  The  kOO   cycle  os- 
cillator is  a  Colpitts  type  with  output  fed  to  the  two  channels  after 
being  amplified  and  shifted  in  phase  90°  with  respect  to  each  other  by 
the  two  RC  type  networks.  The  method  is  similar  to  that  previously  men- 
tioned in  connection  with  the  126  kc  oscillator  and  the  shift  is  in  keeping 
with  the  previous  phase  shifts.  The  purpose  is  to  obtain  the  sura  of  the 
two  channels  across  a  transformer.  It  should  be  mentioned  that  the  shift 
to  U00  cycles  is  to  enable  a  power  gain  to  be  obtained  before  taking  the 
sum  of  the  two  channels. 

Figure  II -6  U00  Cycle  Chopper  and  Amplification  Stages.  -  The  U00  cycle 
chopper  is  a  Model   manufactured  by  the  Bristol  Co.  of  Waterbury,  Conn.. 
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The  necessary  power  amplification  is  provided  by  a  single  stape  in  each 
channel  with  the  two  channels  being  coupled  together  using  a  push-pull 
amplifier.  The  output  to  the  detector  is  the  sum  of  the  two  channels. 
Figure  II-7  Detector  and  Recording  Circuit  -  The  detector  is  a  peak 
reading  type  with  a  time  constant  proportional  to  the  U00  cycle  fre- 
quency. The  Recorder  is  an  Esterline  Angus,  Model  AW,  having  a  1  ma 
movement. 
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3.  OPERATION 

The  technique  of  taking  data  is  quite  simple. 

Previous  calibration  of  the  meter  and  circuits  should  be  made  in 
the  laboratory  using  a  General  Radio  Signal  Generator  805C,  fixed  at  126  kc 
with  zero  modulation.  To  simulate  an  antenna  load  a  value  of  capacitance 
equal  to  the  antenna  capacitance  of  the  particular  aircraft  model  should 
be  shunted  across  the  input  terminals  so  that  the  tuning  capacitor  can  be 
set  for  the  proper  value. 

For  an  input  of  ten  microvolts  the  meter  should  read  a  zero  value 
while  for  a  full  scale  reading  a  value  of  100,000  microvolts  can  be  used. 
Intermediate  points  can  be  located  to  suit  the  user. 

Flights  must  be  at  least  four  hours  long  in  order  to  accumulate  suf- 
ficient data  per  flight.  It  has  been  found  that  a  flight  of  this  duration 
is  necessary  mainly  to  show  the  magnitude  of  the  error  of  the  navigational 
device  for  a  known  course  and  speed  over  a  closed  geographical  loop  of  known 
dimensions. 

Flights  in  the  vicinity  of  thunder  clouds  are  desired  with  a  record  of 
the  time  spent  in  such  locality  marked  on  the  recorder  paper  along  side  of 
the  measurement  record. 
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CHAPTER  III 
RESULTS  AND  CONCLUSIONS 

1.     Results: 

System  II  was  calibrated  in  the  laboratory  and  will  be  delivered  to 
the  Air  Force  some  time  in  April  195U.  No  particular  recommendations  or 
conclusions  can  be  offered  at  this  time  since  no  flights  have  as  yet  been 
made,  but  as  is  typical  with  statistical  data,  the  majority  of  time  spent 
in  investigation  is  used  to  gather  data. 

System  I  was  constructed  and  calibrated  in  the  laboratory  during  the 
month  of  April  195h •  A  unique  testing  set  up,  developed  by  Tanner   was  used 
to  generate  corona  pulses  and  is  illustrated  in  Figure  $. 

Preliminary  recordings  of  the  output  of  this  test  setup  indicated  only 
minute  amounts  of  power  were  being  generated.  It  was  decided  to  take  the 
unit  on  a  test  flight  to  ascertain  whether  any  general  indications  of  pre- 
cipitation static  could  be  obtained. 

A  flight  was  made  on  April  Hi,  195U  from  0830  to  1230  PST  with  route 
as  follows: 

From  Monterey,  California  to  Fresno,  California,  and  then  north-east 
to  the  Sierra  Mountain  chain,  following  this  due  north  to  Yosemite  National 
Park. 

A  long  line  of  cumulus  type  clouds  hung  over  the  ridges  offering  good 
precipitation  static  conditions,  so  the  system  was  placed  in  operation  and 
immediately  the  recorder  indicated  moderate  corona  because  all  discharge 
wicks  were  removed  from  the  plane  previous  to  take-off.  Figures  III  and 
IV  are  parts  of  the  recordings  taken  during  this  flight  and  show  vividly 
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the  difference  in  the  average  power  levels  for  an  area  of  negligible  pre- 
cipitation static  and  an  area  of  heavy  precipitation  static. 

Conclusions: 

Investigations  of  precipitation  static  made  using  System  II  will 
record  values  related  only  to  the  average  power  input  to  a  particular 
device  so  it  is  obvious  that  this  information  is  of  no  value  to  anyone 
interested  in  the  actual  power  spectrum  or  the  possible  frequency  distri- 
butions of  precipitation  static.  It  was  for  this  purpose  that  System  I 
was  designed. 

Since  the  data  gathering  period  will  take  considerable  time,  only 
a  brief  outline  of  the  method  and  approach  will  be  given. 

The  approach  to  t  he  problem  can  be  again  allied  to  the  familiar  treat- 
ment of  fluctuation  noise,  using  the  autocorrelation  function  and  its 
relation  to  the  Fourier  cosine  integral. 

The  autocorrelation  function  /  (T)  of  a  time  function  i(t)  is  defined 
as  the  time  average  of  the  product  of  all  pairs  of  points  i(t)  separated 
by  an  interval  .  .  *T  , 

Substituting  from  Equation  1-20  - 

p(rO~  v  [°° c(t)  c  (t-hr)<Ltt  Tit)1" 

In  the  above  the  second  term  is  the  square  of  the  d-c  value  I  (t). 

If  a  current  i(t)   is  defined  as  the  fluctuation  of  the  current  about  its 

mean  value,  then 
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i  (t)  =  l(t)  -  I 

^(T)  =  v  \   i(t)   i  (t  +  T) 

^(0)  s  i  ^  =  v  V  i  (t)  dt 
Note  that  the  total  noise  power  depends  not  only  upon  the  average  current 
but  on  the  shape  of  the  current  pulse  as  well. 

The  noise  power  spectrum  G.  (w)  can  be  obtained  from  the  autocorre- 
lation function  by  the  Fourier  cosine  integral  relation,  and  is  found  to 
be 

Gi(w)  =  2  v     k   e     (w  To.)2   2(1  -  cos  w7„-  w  sin  w  TJ 
(  w  T+)k 

where  e  is  equal  to  the  charpe  of  electron,  and  can  be  plotted  as  shown 
in  Figure  6. 

The  noise  power  in  the  frequency  range  f,  to  f«  is 
\  G(f)  df 

In  the  frequency  range  for  which  wTt^^  ,  G(f)  is  nearly  constant 


in  value.  This  is  true  no  matter  what  the  shape  of  the  current  pulse  is 
since  the  value  of  the  Fourier  Integral  at  zero  frequency  depends  only 
upon  the  zero  moment  of  the  function.  Thus  G(0)  depends  only  upon  the 
area  under  f^   (t)  and  not  upon  its  shape,  so  that  the  noise  power  in  a 
bandwidth  L  f    is2  ve^f. 
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FIGURE  6 
Plot  of  Noise  Power  Spectrum  versus  uj'T 
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A  device  then,  having  a  bandwidth   that  is  less  tnan  one  radian  can  di- 
rectly measure  increments  of  the  power  spectrum.  This  work  is  now  in  pro- 
gress using  System  I  using  the  noise  model  described  in  the  beginning  of 
this  chapter.  It  is  obvious  that  the  noise  model  device  is  inadequate  as 
a  generator  for  heavy  corona  conditions  and  investigation  and  design  of  a 
superior  device  is  contemplated. 
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APPENDIX  I 
THE  NORMAL  DISTRIBUTION 
1.  Probability: 

(a)  If  an  event  is  certain  to  occur  it  is  assigned  an  ar- 
bitrary value  of  unity.  For  example: 

In  the  toss  of  a  single  die  the  probability  of 


a  six    is   1/6 

a  five 

•   1/6 

a  four 

'   1/6 

a  three 

'   1/6 

a  two 

•   1/6 

a  one 

•   1/6 

thus,  if  six  probabilities  are  mutually  exclusive,  the  probability  that 
either  one  or  the  others  will  occur  is  equal  to  the  sum  of  the  separate 
probabilities. 

(b)  When  a  particular  result  in  any  trial  is  independent  of 
the  results  obtained  in  any  other  trials,  is  called  a  random  process.  For 
instance: 

A  large  number  of  independent  trials  in  tossing  a  coin  might  reveal  a 
consecutive  string  of  heads.  This  will  not  affect  the  next  toss  since  the 
probability  is  still  1/2  that  the  coin  will  come  up  either  a  head  or  a  tail. 

(c)  The  probability  of  the  simultaneous  occurrence  of  n  events 
is  the  product  of  the  n  individual  probabilities  of  n  events.  Thus  the 
probability  of  winning  all  tosses  in  a  season  of  eight  footballs  games 

o 

would  be  (0.5)     =     0.0039 

The  above  three  rules  form  the  basis  of  the  probability  distribution 
laws  upon  which  the  mathematical  description  of  noise  is  based. 
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2.  Distribution  formulas: 

If  the  symbol  P_(n)   is  made  to  stand  far  the  probability  that  an 
event  will  occur  exactly  n  times  in  m  independent  trials,  then  the  basic 
distribution  formula,   (Bernoulli  Distribution),   is 


t,  /    \     -     J"1      n,  vm-n  n 

Pm(n)     "     Cn    p  (1  -  p)         =  ml  pn  (1  -  p) 

n  I   (m  -  nj    I 


m-n 


where  p  is  the  probability  for  the  event  to  occur  in  any  particular  trial. 
This  specializes  to,   (Poisson  Distribution), 


* 


p  /    \  -        -         /-Nn  for      m   »        n 

Pm(n)  -     e-n         (n  p   ^       1 

ni  m    >>       1 


where 


n     =    mp 

When  n  becomes  very  large  as  in  the  case  of  high  density  Poisson 
v    noise,  the  above  further  specializes  to  the  Gaussian  Distribution. 


Pt>>     =  1  e  -(n-H>2/  2* 

when  n    >>   1 


n    -     n      ^<  1 


n 
n  -  H      >>      i 
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3»  The  Central  Limit  Theorem 

The  Central  Limit  Theorem  states  that  subject  to  very  lenient 
conditions,  and  regardless  of  the  amplitude  distributions  of  the  indi- 
vidual components,  as  the  number  of  components  n  approaches  infinity, 
the  amplitude  distribution  of  the  current  approaches  the  Gaussian  or 
normal  distribution. 

Thus  it  is  possible  to  relate  the  limiting  form  of  the  Gaussian 
Distribution  to  the  Mathematical  Description  of  Noise  (Chapterl)  by 
virtue  of  the  large  number  of  independent  events  (electron  crossings) 
which  add  up  to  give  the  independent  variable  I.  This  distribution 
will  be  independent  of  the  amplitude  distributions  of  the  individual 
events  (shapes  of  the  electron  pulses).  It  is  required,  generally  speak- 
ing, that  the  contributions  of  the  individual  components  be  small  compared 
with  the  total  value  of  the  sum 
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APPENDIX  II 
PROPERTIES  OF  NARROW  BAND  NOISE 

1.  Examining  the  statistical  properties  of  the  output  of  a  narrow 
band  filter  plotted  in  Figure  A -I,  and  the  assumption,  that  the  band 
width  of  the  filter,   f^^  To  >the  output  can  be  expressed  in  the  form 

e(t)  *  E(t)     cos     (wmt-*V(t))  II-l 

provided  that  the  input  noise  has  a  Gaussian  distribution. 

E(t)  then  represents  the  magnitude  of  the  envelope  of  the  output  and 

^(t)  the  phases  changes. 

Using  the  Fourier  representation  of  the  output  as 

e(t)»y  (a  cos  nw  t  +  b  sin  nw  t)  II-2 

<—      n      o     n       o 


where  the  summation  is  to  include  all  frequencies  which  is  in  the  range 
of  the  filter. 

Setting  up  the  sum  and  difference  terms  so  as  to  show  relationship 
to  a  sine  wave  of  frequency  w  :• 

cos  nwQt  =  cos  T  (nwQ  -wm)t  +■  w  t  \  II-3 

=    cos       (nw    -w  )t  cos  w  t  -  sin  (nw    -w  )  sin  w  t 
o       ra  m  o       m  m 

sin  nw  t     =     sin  I   (nw     -w  )t+wt  II-U 

v  \         o       m  m     I 

=     sin       (nw      -    w  )t  cos  w  t  4-  cos  (nwrt  -  w„)t  sin  w  t 


& 


e(t) 


m 
Frequency 

FIGURE  A- I 
Output  of  a  Narrow  Band  Filter 


£6 


Letting   e    (t)  represent  the  summation  of  all  the  coefficients  of 

cos  wmt  in  the  expression  for  e(t),  and  egin(t)  represent  the  summation 

of  the  coefficients  of  sin  w  t,  then 

m  ' 


(t)  S  I  fcn  cos  (nwo  "  wm>t  +  bn  sin  ^nw0  "  *»>*]    ll* 


e 

cos 


and 


eein<t}  =  l_  Lbn  C0S  (nWo  -  V*  -  an  sin  ^o  "  V*  J    U"6 


v\ 


Substituting 

e(t)  =     ecos  (t)   cos  V  +   esin  (t)   Sin  V  IL"7 

Thus,  e         (t)     and  e   .      (t),  are  slowly  varying  compared  with 
cos  wm  t  or  sin  wmt. 


From  the  Expression  III-l,  the  magnitude 


E(t)  = 


eL     U)  +  «?,-     (t)  H-8 


^  cos     N  sin 

and 

rf(t)«     tan"1     -esin  (t)  /  ecog  (t)  II-9 

The  random  variables  ec  (t)  and  e  .  (t)  are  sums  of  independent 
variables  with  Gaussian  amplitude  distributions.  The  average  power  can 
then  be  written  as 
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E 


e  e    .  2  11-10 

cos  sin  e 


where 

2 


and 


e 
cos 


=  1    4 


n 


e2.       =Z.         b2 

sin  n  n 
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